The prospects for a new search for the electron electric dipole moment in solid 
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We address a number of issues regarding solid state electron electric dipole moment (EDM) 
experiments, focusing on gadolinium iron garnet (abbreviated GdIG, chemical formula Gd3Fe50i2) 
as a possible sample material. GdIG maintains its high magnetic susceptibility down to 4.2 K, which 
enhances the EDM-induced magnetization of a sample placed in an electric field. We estimate that 
lattice polarizability gives rise to an EDM enhancement factor of approximately 20. We also calculate 
the effect of the demagnetizing field for various sample geometries and permeabilities. Measurements 
of intrinsic GdIG magnetization noise are presented, and the fiuctuation-dissipation theorem is 
used to compare our data with the measurements of the imaginary part of GdIG permeability 
at 4.2 K, showing good agreement above frequencies of a few hertz. We also observe how the 
demagnetizing field suppresses the noise-induced magnetic fiux, confirming our calculations. The 
statistical sensitivity of an EDM search based on a solid GdIG sample is estimated to be on the 
same level as the present experimental limit. Such a measurement would be valuable, given the 
completely different methods and systematics involved. The most significant systematics in such an 
experiment are the magnetic hysteresis and the magneto-electric effect. Our analysis shows that it 
should be possible to control these at the level of statistical sensitivity. 



O 

• ^ ■ 

^ ; 

Oh: 



> 
d 

OO 

o 



I. INTRODUCTION 

One of the ways to test the Standard Model in low- 
energy experiments is to search for parity and time- 
reversal- violating permanent electric dipole moments of 
particles such as the electron and the neutron. It is 
well-known that the observed matter-antimatter asym- 
metry of the Universe can not be explained by the CP- 
violation mechanisms within the Standard Model 
A variety of theories suggesting new sources of CP- 
violation exist, many of them are already constrained 
by the current experimental limit on the electron EDM: 
dp < 1.6 X 10~^^ e-cm A number of experimental 
EDM searches are currently under way, or being devel- 
oped, they involve diatomic molecules 0, 0, @, molec- 
ular ions al, cold atoms f?\, liquids and solid-state 
systems [9|. 

Following the original suggestion of Shapiro [l^l, we 
have been investigating the possibility of an improved 
limit on the permanent electric dipole moment (EDM) 
of the electron by use of paramagnetic insulating solids 
together with modern magnetometry [ll| . The basic idea 
is that when a paramagnetic insulating solid sample is 
subjected to an electric field, if the constituent atoms or 
ions of the solid have an electric dipole moment, the atom 
or ions spins will tend to become spin-polarized. Because 
the atoms or ions carry a magnetic moment, the sample 
acquires a net magnetization, which can be detected by 
a SQUID magnetometer inside a set of superconducting 
magnetic shields in a liquid helium bath (see Fig.[T]). 

Three solid state electron EDM experiments have been 
reported in literature. The first of these used a nickel- 
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FIG. 1: A schematic experimental setup for an EDM search. 



zinc ferrite, obtaining a limit of about 10~^^ e-cm 
Another experiment, described in Ref. [l^, obtained a 
limit that is lower by a factor of 40. This experiment 
employed a converse effect, in which a sample electric 
polarization appears when the sample is magnetized. 
The third experiment employs Gadolinium Gallium Gar- 
net (GdGG), which has previously been identified as a 
promis ing paramagnetic system for this type of exper- 
iment jl4| . We have performed preliminary measure- 
ments, indicating that the anticipated sensitivity can be 
achieved, pending a 1/T behavior of the paramagnetic 
susceptibility to temperatures of order 10 mK. However, 
attaining an improved electron EDM limit requires that 
GdGG remains a simple paramagnetic system and that 
spin-glass or other effects do not enter down to tempera- 
tures of order 10 mK. Indeed, it is now known that such 
effects become very important near 0.2 K [l^, where the 
magnetic susceptibility reaches a peak value of x = 0.15 
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(CGS units). This value is actually larger than the sus- 
ceptibility given by the Curie-Weiss law, x « {T — Tc)~^, 
which describes the susceptibility of GdGG above 1 K, 
where Tc ~ —2 K for GdGG. We are continuing to inves- 
tigate the ultimate sensitivity of a GdGG-based exper- 
iment, and the relaxation time of the spins remains to 
be determined. It is possible that this time will be long 
enough that extra noise due to spin fluctuations will be a 
limitation jl6j . as we discuss below. In the meantime we 
have mounted a search for other materials whose proper- 
ties make them attractive for conducting this type of an 
EDM experiment. 



II. THE SEARCH FOR THE BEST MATERIAL 

A. The magnetic susceptibility 

For a solid state electron EDM experiment, materials 
with large magnetic susceptibility (or relative permeabil- 
ity) are desired. This is because, neglecting boundary ef- 
fects (see Section IIII| , the EDM- induced magnetization 
is given by 

M = xkdeE/na, (1) 

where x is the magnetic susceptibility, k is the effective 
EDM enhancement factor in the material under study, 
is the electron EDM, E is the applied electric field, and 
jia is the magnetic moment of the paramagnetic atoms or 
ions that carry the electron EDM {^a = 8^3 for Gd+'', 
/is is the Bohr magneton). This assumes that all of the 
magnetization is due to the magnetic moments of the 
heavy atoms, that carry an EDM (a good approximation 
for all cases considered here). We see immediately that 
X should be as large as possible. Note that the sample 
temperature does not explicitly enter into this formula. 

Garnet ferrites with rare earth substitution can have 
large permeabilities at temperatures below 100 K. We 
have studied the complex permeability (/i = /i' — i/i", 
/i = 47rx) of mixed Gadolinium- Yittrium iron garnets 
(GdYIG, chemical formula Gd(3_2;)Ya;Fe50i2) down to 
2 K . Our data show that for pure GdIG ceramic, fi' — 
77 at 4.2 K, and for mixed Gdi.8Yi.2Fe50i2 ceramic, fi' — 
50 at 4.2 K. The immediate conclusion is that at 4.2 K, 
in the absence of sample-shape effects, the EDM-induced 
magnetization in pure GdIG is enhanced by a factor of 
more than 300 compared to that in GdGG, which has x = 
0.016 at 4.2 K [I3. The present manuscript describes 
our study of the feasibility of an EDM experiment with 
GdIG. 



B. The effective EDM enhancement factor and the 
internal electric field 

The enhancement factor k that appears in Eq. ^ de- 
termines the scale of the effective electric field acting on 



the Gd+3 ion EDM in GdIG. To calculate this, let us con- 
sider the energy shift S resulting from a non-zero EDM 
of an ion in the lattice of an ionic solid. A naive estimate 
of this energy shift can be written down as follows: 

e + 2 

Snaivc ~ —diE\oca\ ~ — diE, (2) 

where di is the Gd^^^ ion EDM, i?iocai is the local electric 
field at a Gd+'^ ion site, e is the dielectric constant of 
GdIG, and E is the applied electric field. Here we used 
the Lorentz relation for the local field in a dielectric (in 
the approximation of a cubic crystal structure) [18]. The 
magnitude of dj induced by the electron EDM is calcu- 
lated in Ref. [l||: Id,] « 2.2d^^ . With e = 15, we get the 
naive estimate (5naivc ~ —Hd^E. 

A more rigorous calculation is presented in Refs. [H, 
|20| . [2]| . The energy shift is expressed in terms of the 
displacement x of the EDM-carrying ion with respect to 
its equilibrium position in the unit cell: 

(5 w -0.1 — — • 27.2 eV, (3) 

as eos 

where as is the Bohr radius, e is electron charge, and 
27.2 eV is the atomic unit of energy. The Gd+'^ ion dis- 
placement X is related to the applied electric field through 
the dielectric constant: x — {eion ~ l)E/12TTen, where 
£ion is the dielectric constant contribution due to ion dis- 
placement and n « 10^^ cm~'^ is the Gd density. We are 
assuming that the displacement of Fe ions is much smaller 
due to the stiffer Fe-0 bond. The other contribution to 
the dielectric constant in GdIG is the "electronic polar- 
izability" - the polarization of the individual ions in the 
electric field. We can estimate this contribution from the 
refractive index n w 2.2 of GdIG at optical frequencies, 
where the lattice polarizability does not contribute. Thus 
£ion ~ 15 — 2.2^ « 10 is a reasonable estimate for the lat- 
tice polarizability. With this value, the EDM-induced 
energy shift in GdIG is: 

5 w ~2MeE. (4) 

Thus the effective enhancement factor is fc = 20. As 
pointed out in Ref. [l^ . this is a conservative estimate, 
and the true result may be up to a factor of two larger. 
Note that the final value of k is quite close to the naive 
estimate discussed above. 



III. THE EFFECT OF THE SAMPLE SHAPE 

Above we described two enhancement mechanisms for 
an EDM-induced signal in GdYIG ferrites: large mag- 
netic permeability and local electric field enhancement 



This naive estimate of the enhancement factor gives the incorrect 
sign if the actual value di —2.2de from Ref. [l9j is used. The 
issue of this sign is addressed in Ref. |20( 
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by the crystal lattice. The combination of these makes 
this an attractive system for mounting an EDM search. 
It must be noted, however, that sample geometry- 
dependent demagnetizing fields eliminate part of the ad- 
vantage of having a large-permeability material. It is 
well known, for example, that it is difficult to magne- 
tize a thin disk of large-permeability magnetic material 
along its axis of symmetry (permanent magnets with such 
magnetization exist because their magnetization is satu- 
rated, i.e. their linear permeability is close to 1). This 
is explained by the Maxwell's equation V • B = 0, com- 
bined with the relation B ^ H + 4ttM (CGS units). If 
a sample has some macroscopic magnetization M, then 
on its surface V ■ H — — 47rV • M ^ gives rise to 
the "demagnetizing field" if, which opposes the magne- 
tization M . The magnitude of this demagnetizing field 
depends on the sample geometry and permeability. For 
a high-permeability material the resultant magnetization 
can be much less than that in the case of no demagnetiz- 
ing field, which happens for a long thin needle geometry, 
for example. 

The demagnetizing field, and the sample magnetiza- 
tion, can be calculated exactly when the material is in 
the form of an ellipsoid [l^l, but not for a disk-shaped 
sample, which would be used an an EDM search (some 
approximate expressions for saturated hard ferromagnets 
are obtained in Refs. l23i, 12^, |25|, i26,] ) . We used the finite- 
difference method to numerically solve Maxwell's equa- 
tions V -B — Q and V x Jf = 0, with a cylindrical sample 
of radius r, thickness h, and permeability /i. In such a 
geometry the magnetization is not uniform, but we can 
calculate the total magnetic flux $ through a circular 
loop positioned on the surface of the cylinder halfway 
between its top and bottom. When /i ^ r, there is no 
demagnetizing field, and <& — AnMnr'^, where M is the 
EDM-induced magnetization given by Eq. ([!]). In gen- 
eral, however, we have to write 

$ = 5 • 47rAf7rr^, (5) 

where g — g{h/r,ii) is the flux suppression geometric 
factor that depends of the geometry and the permeability 
of the sample. The results of our numerical calculations, 
expressed as g{h/r,fj,), are shown in Fig.[2l 

There is clearly a significant suppression of the EDM- 
induced magnetic flux for a reasonable sample geometry 
(h/r ~ 1) and for ^ ~ 50. Certainly the suppression 
factor g does not give the whole story, since there are 
several other geometry-dependent factors that enter into 
the final EDM sensitivity, such as the pickup loop area 
and inductance, as well as the electric field for a given 
experimentally achievable voltage applied to the sample. 
Further discussion of these factors can be found in Sec- 
tion |Vl It is also important to note that for h/r ^ 1, 
there is a significant suppression even for — 1 <C 1, as 
would be the case in a search for a nuclear Schiff moment, 
for example Q. Therefore, in any magnetization-type 
EDM search, a thin disk-shaped sample would lead to a 
loss in sensitivity, instead samples with /i ~ r must be 
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FIG. 2: The magnetic flux suppression geometric factor for a 
permeable cylindrical sample. 

used^. 



IV. THE MAGNETIC NOISE MEASUREMENTS 

Even with a substantial signal suppression by the de- 
magnetizing fields, it is possible, by carefully choosing 
the sample dimensions, to design an EDM search of 
formidable sensitivity. Before a full-scale experiment 
with Gd-containing ferrites is attempted, however, the 
intrinsic magnetic noise of the material must be mea- 
sured. If this noise is greater than the magnetometer 
noise, the sensitivity of the EDM-induced magnetization 
measurement is reduced [l6j . 

The setup used for the GdIG and GdYIG magnetic 
noise measurements is shown in Fig. [31 The sample is 
a sintered ceramic toroid with a square cross-section, 
outer diameter 3 cm, inner diameter 1 cm, and thickness 
1 cm. Experiments were performed with ceramic sam- 
ples of Gd3Fe50i2, obtained from Pacific Ceramics ^2^. 
The superconducting pickup loop was made from enamel- 
insulated niobium wire of thickness 0.003 in. It was con- 
nected to the input terminals of a Quantum Design DC 
thin-film SQUID sensor (model 50). The SQUIDs were 
connected to the control unit (model 5000) via 4-meter 
MicroPREAMP cables. The measurement bandwidth 
was limited to 1 kHz by an analog filter in the SQUID 
controller, to prevent aliasing. The intrinsic flux noise of 
the SQUID magnetometers in this frequency range was 



^ This conclusion is also valid if an atomic magnetometer (or some 
other magnetometer) is used to detect the magnetic field outside 
the sample, instead of a SQUID that detects the magnetic flux 
through a pickup loop. 
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2 /i$o / ■\/Hz, which corresponds to magnetic field noise of 
several fTZ-s/Hz, the exact value depending on the area 
of each pickup loop. This noise level was always much 
less than the measured sample magnetization noise. The 
SQUID noise 1/f corner was around 0.1 Hz. The SQUIDs 
and the sample were mounted on G-10 holders, carefully 
fixed to minimize vibrations, and enclosed in supercon- 
ducting shielding, made of 0.1-mm thick Pb foil, glued to 
the inner surfaces of two G-10 cylinders. The magnetic 
field shielding factor of our setup was greater than 10^. 
The entire setup was submerged in liquid helium inside a 
Janis model lOCNDT dewar. To dampen out vibrations, 
the dewar was mounted on rubber supports and placed 
on a 15" diameter support ring, comprising a partially in- 
flated bicycle inner tube; however the residual vibrations 
of the SQUID pickup loops with the respect to the resid- 
ual magnetic field gradients'^ inside the superconducting 



magnetic shields were enough to cause the pronounced 
high-frequency (between 10 Hz and 1 kHz) peaks in the 
SQUID measurements shown in Fig.[5l Temperature was 
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FIG. 3: The experimental setup for magnetic noise measure- 
ments. 



FIG. 4: The three SQUID pickup loop geometries used for the 
magnetic noise measurements: a. pickup loop in the shape of 
figure-eight for external magnetic field rejection; b. U-shaped 
pickup loop for rejection of magnetic field lines internal to 
the toroid; c. pickup loop wound around the toroid perime- 
ter to detect axial magnetic field. Magnetic field vectors are 
schematically shown in red. Bo is the magnetic field outside 
the toroid. Bi labels an azimuthal magnetic field line, looping 
around inside the toroid. 



^ These gradients are probably caused by the magnetic flux 



trapped at some pinholes in the Pb foil. 
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measured with a LakeShore model DT-670C-SD silicon 
diode temperatm'e sensor, accm'ate to 0.1 K. 

Magnetic noise measurements were done with several 
pickup-loop geometries to investigate the effect of the 
demagnetizing field on the noise. It is energetically fa- 
vorable for the magnetic field lines to stay inside a high- 
permeability material, therefore most of the flux due to 
the magnetic noise in the material is confined inside the 
toroid, as indicated in Fig. [D This flux can be detected 
by a SQUID pickup loop wound in the form of a figure- 
eight, i.e. in the gradiometer configuration (Fig.[3|i). The 
spectrum of the measured magnetic field noise for this ge- 
ometry is shown as the red curve marked (a) in Fig.O In 
this case, the magnetic field lines were closed inside the 
toroid, and, to a good approximation, there was no de- 
magnetizing field (geometric factor g = 1), therefore this 
was a measurement of the intrinsic magnetization noise 
of the ferrite sample. 

To demonstrate that the observed magnetic field noise 
is due to the magnetization inside the material, rather 
than the magnetic field noise originating outside the sam- 
ple, we performed noise measurements with a SQUID 
pickup loop wound in the U-shaped form, i.e. in the 
magnetometer configuration perpendicular to the plane 
of the toroid (Fig. |3|d). The spectrum of the measured 
magnetic field noise for this geometry is shown as the 
green curve marked (b) in Fig. [S] In this case, the flux 
due to the magnetic field lines that were closed inside the 
toroid was not coupled to the SQUID, and only the mag- 
netic field that entered or exited the toroid was detected. 

Finally, to simulate the geometry of an EDM exper- 
iment as closely as possible, we measured the magnetic 
field noise that couples to the SQUID with the pickup 
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FIG. 5: The results of the magnetic noise measurements. The 
curves are marked with the corresponding pickup loop geome- 
try, shown in Fig. U) Two noise expected from the fluctuation- 
dissipation theorem (FDT) for geometries (a) and (c) is also 
shown. 



loop wound around the perimeter of the sample, as shown 
in Fig.[4j:. In this case, the magnetic fiux noise was given 
by the GdIG magnetization noise, reduced by the ap- 
propriate geometric factor that takes into account the 
demagnetizing field, as in Eq. ([5]). The results arc shown 
as the blue curve marked (c) in Fig. [S] The pronounced 
peaks in all three of the noise data sets are caused by 
vibrations, as discussed above. 

Let us proceed to the analysis and discussion of these 
magnetic noise data. The magnetization noise in a mate- 
rial with a known complex permeability ^ = fj! — i^" can 
be calculated using the fluctuation-dissipation theorem 
(FDT) [13: 

(M2)^ = -i-M:^"(^), (6) 
ZnV LO 

where {M'^)^^ is the power spectral density of the mag- 
netization noise at frequency uj, V is the sample vol- 
ume, ks is Boltzmann's constant, and T is temperature. 
Complex permeability has been measured in GdIG and 
mixed GdYIG ceramics in a range of temperatures (2 K 
to 295 K) and frequencies (100 Hz to 200 MHz) p?]. In 
order to compare our noise measurements with the FDT, 
we measured the complex permeability of GdIG at 4.2 K 
in the frequency range of 0.01 Hz to 1 kHz. This was 
done by winding two 23-turn excitation coils in series, one 
around the GdIG toroid, and one around a toroid of the 
same dimensions, but made of G-10. Each of the toroids 
was connected to a SQUID via a one-turn pickup loop. 
The GdIG complex permeability was extracted from the 
relative amplitude and phase shift of the SQUID signals, 
as a sinusoidal current with frequency ranging between 
0.01 Hz and 1 kHz was applied to the excitation coil. 
The resulting values of /i" were used in Eq. © to gen- 
erate the magnetic field noise expected from the FDT. 
This is shown by circles in Fig. [5l in the case of no de- 
magnetizing fields (geometric factor g = 1), which is a 
good approximating for geometry (a) . 

For frequencies greater than approximately 5 Hz, there 
is very good agreement with the magnetic field noise mea- 
sured with the SQUID for this geometry. We interpret 
this as a verification of the validity of the FDT for our 
system in this frequency range. For lower frequencies, 
however, the FDT seems to systematically underestimate 
the magnetic field noise, by as much as two orders of mag- 
nitude. We speculate that the excess magnetic field noise 
is due to very slowly-relaxing degrees of freedom, which 
come to thermal equilibrium on the time scales much 
longer than the duration of our experimental run (a few 
hours) [29} . For such degrees of freedom, the effective 
noise temperature that should be used in Eq. ([6]) is much 
greater than 4.2 K, thus they contribute greater noise. 

In Section [m] we discussed how the sample shape 
affects the EDM signal, presumably the magnetization 
noise is also reduced by the same demagnetizing fields. 
This is confirmed by the measurements of the noise with 
the SQUID pickup loop around the sample perimeter 
(Fig St)- Indeed, if the FDT prediction © is muhi- 
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plied by the geometric factor appropriate for the sample 
dimensions {h/r — 0.67 resulting g — 0.017)^, the re- 
sult is the expected magnetic noise in the geometry of 
FigSh- This is shown by the squares in FigO and once 
again, the agreement with the SQUID measurements is 
remarkable above frequencies on the order of a few Hz. 
We can therefore conclude that we predict the magnetic 
field noise in any geometry by taking into account the 
geometric factor, as described in Section Hill 

V. THE PROJECTED SENSITIVITY OF AN 
EXPERIMENT WITH GDIG 

Having measured the magnetization noise in GdIG, we 
can estimate the statistical sensitivity of an EDM exper- 
iment with a GdIG sample. Combining Eqs. ([T]) and (O, 
we obtain the EDM-induced magnetic flux through a 
pickup loop around a cylindrical sample: 

2 . kd^E . . 

$EDM = ffTTr • 47rx-3 , (7) 

where all the symbols have been defined previously. The 
statistical sensitivity is limited by the material magne- 
tization noise, given by Eq. ([S]), the resulting magnetic 
flux noise through the pickup loop is given by: 

2 I 1 fesT/i" 

where / = u!/2tt is the electric field reversal frequency. 
We can now express the statistical sensitivity in terms of 
experimental parameters: 

. , / h ksTn" 

^'^"^ — t:tf\~ — ~f — ' 

TTxkVo y 7rr^ / 

where we introduced the amplitude Vq of the voltage ap- 
plied to the electrodes. Interestingly, the geometric factor 
g does not enter this sensitivity estimate, this is because 
it appears in both the EDM-induced flux and the noise- 
induced flux. If the magnetization noise had turned out 
to be lower than the intrinsic noise of the SQUID, the geo- 
metric factor would have entered into the sensitivity esti- 
mate (so would the SQUID sensitivity and the number of 
turns of the SQUID pickup loop, see Appendix [X)) . Equa- 
tion @ assumes a disc-shaped geometry, but the sample 
dimensions have not been specifled. It is clear that the 
best sample is a thin disc of large radius: decreasing the 
thickness h increases the applied electric field for a given 
voltage Vq , and increasing the radius r increases the sam- 
ple volume, which decreases the magnetization noise, ac- 
cording to Eq. It is also clearly advantageous to 



* We approximate the sample as a disk. Our simulations show that 
the presence of the central hole does not change the geometric 
factor by more than a few percent 



reverse the electric field as fast as possible, however dis- 
placement currents limit the maximum realistic reversal 
frequency to 10 Hz (see below for more detailed discus- 
sion of systematics) . For a quantitative estimate of the 
EDM sensitivity we take the following experimental pa- 
rameters: Vf) — 10 kV, h — 0.1 cm, r — 3 cm, / = 10 Hz, 
resulting in 5de ~ 4 x 10"^"' e-cm/VHz. After 10 days of 
integration, a statistical sensitivity at the level of 

(54 =i 3 X 10"^^ e-cm (10) 

can be achieved. 

A number of systematic concerns for a solid state 
sample-based EDM search have been outlined in 
Ref. [ll|. With a ferromagnetic system, magnetic hys- 
teresis is also crucially important. The main concern is 
electric field-correlated sample magnetization caused by 
the magnetic fields generated by displacement currents. 
It should be noted that these magnetic fields tend to be 
perpendicular to the EDM-induced magnetization, but 
some degree of misalignment is inevitable. The current 
associated with applying 10 kV to a 300 pF sample at 
10 Hz reversal frequency is about 30 /lA, which gives rise 
to magnetic field of order — 2 x 10^^ A/cm. The 
resulting remanent magnetization Md, present after the 
current is brought back down to zero, can be calculated 
from the parameters of the Raleigh hysteresis loops mea- 
sured for GdIG in Ref. [l3| • Taking a reduction of a factor 
of 10 due to perpendicularity, the resulting magnetization 
along the electric field is Aljj ~ 10^^^ emu/cm"^. Accord- 
ing to Eq. (I5|), the induced magnetic flux is suppressed 
by the geometric factor, which, for a disc of ft- = 0.1 cm 
and r = 3 cm, is approximately 10~^. The resulting sys- 
tematic is at the level of statistical sensitivity given in 

Eq. cni). 

Another possible systematic is the so-called magneto- 
electric effect [l^j. The applied electric fleld induces a 
strain in the sample (electrostriction), which then in- 
duces a magnetization correlated with the electric field, 
provided the sample already possesses some non-zero 
magnetization (inverse magnetostriction). This effect is 
described by the term ^E'^H^ (we omit the tensor in- 
dices) in the electromagnetic free energy, which is al- 
lowed by both the parity and the time-reversal symme- 
tries. This has been observed in YIG ^ and GdIG 
The "magneto-electric susceptibility" Pm is defined in 
the following manner: M — /3m{HQ)E'^ , where Hq is 
the external magnetic field, E is the applied electric 
field, and M is the resulting magnetization. For GdIG 
single crystals in the low-field limit {Hq < 30 Oe), 
the magneto-electric susceptibility was measured to be: 
l3m{Ho) — 10^^ • Hni^e) emu/cm'^, where -ffo is ex- 
pressed in Oestreds [31j. For a poly-crystalline sample, 
this should be reduced by a factor of 3 (the average 
over crystallite orientations). With the electric field of 
100 kV/cm, and assuming 1% reversal accuracy, the ex- 
ternal magnetic field Hq has to be kept below 3 x 10~^ Oe 
in order to control this systematic at the level of statisti- 
cal sensitivity. This appears feasible with a combination 
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of Metglas magnetic shielding and superconducting lead 
shields. 

Finally, we point out that, if a gradiometer setup is 
used with a large-permeability sample, the gradiometer 
tuning has to be adjusted, see Appendix IB] for analysis. 



VI. CONCLUSION AND OUTLOOK 

We have performed a detailed feasibility study for an 
electron electric dipole moment search with a solid ce- 
ramic sample of Gadolinium iron garnet. This material 
is attractive because it is an excellent insulator, and it 
maintains a high magnetic permeability down to liquid 
helium temperatures, fj, — 77 at 4.2 K, which enhances 
the EDM-induced magnetization and allows the design of 
an experiment based on SQUID-magnetometer detection 
and superconducting magnetic shielding. In addition, 
the electric field acting on the electron EDM is ampli- 
fied by the crystalline lattice polarizability, resulting in 
an effective EDM enhancement factor of 20. When cal- 
culating the EDM-induced magnetic flux, one must take 
into account the demagnetizing field, which suppresses 
it by a factor that depends on the sample geometry and 
permeability. The dominant sensitivity limitation in the 
proposed experiment, however, is the intrinsic magneti- 
zation noise of the Gadolinium iron garnet samples at 
4.2 K. We measured this noise and verified that our data 
are consistent with the fluctuation-dissipation theorem 
in the 5 Hz to 1 kHz frequency range. At lower fre- 
quencies the noise seems to be dominated by magnetic 
degrees of freedom that are out of equilibrium with the 
liquid helium heat bath. We also verifled the expected 
scaling of the noise-induced magnetic flux with the sam- 
ple geometry. We estimate the statistical sensitivity of 
an EDM search based on a solid GdIG sample to be 
3 X 10^^^ c-cm after 10 days of averaging. This is slightly 
above the present experimental limit, but such a mea- 
surement would still be highly valuable, given the com- 
pletely different methods and systematics involved. We 
also estimate the most likely systematics: the magnetic 
hysteresis and the magneto-electric effect. Our analysis 
shows that it should be possible to control these at the 
level of statistical sensitivity. 

It seems feasible that the magnetic noise can be re- 
duced by diluting the GdIG powder with a non-magnetic 
mixer (such as teflon) prior to pressing the ceramic sam- 
ples, Another possibility is substituting some of the Fe 
ions in the lattice with paramagnetic ions, such as Ga, 
making mixed Gd3Fe(5_^)Gaa;Oi2 ceramics. This will 
probably reduce the magnetic susceptibility as well as the 
magnetization noise, but, by carefully choosing the sam- 
ple dimensions, this loss can be compensated by an im- 
proved geometric factor. Yittrium-substituted Gadolin- 
ium garnets with chemical formula Gd3_2:Y2,Fe50i2 have 
been studied in Ref. [l^, and the measured complex 
permeability at 4.2 K leads to a poorer EDM sensi- 
tivity estimate for such ceramics. We have, in addi- 



tion, identified an extremely promising material where 
the sensitivity gain is provided by the crystal lattice, 
rather than by the large magnetic permeability [s^ . 
This material is gadolinium molybdate (chemical formula 
Gd2(Mo04)3), which is a ferroelectric-ferroelastic with 
Curie temperature of 160°C and spontaneous polariza- 
tion of 0.25 /iC/cm^. Preliminary estimates and mea- 
surements indicate that a two-order of magnitude im- 
provement over the present limit is possible with this 
material. We are currently pursuing further studies of 
Gd2(Mo04)3, as well as other magnetic ferroelectrics. 
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APPENDIX A: THE MULTI-TURN SQUID 
PICKUP LOOP 



If the sample magnetization noise happens to be below 
the SQUID noise, the EDM sensitivity is limited by the 
SQUID noise. To achieve the best coupling of the sample 
flux to the SQUID, a multi-turn pickup loop may be used. 
This appendix shows that having a several-turn pickup 
loop may be advantageous, and suggests how the number 
of turns may be optimized. 

Let us denote by N the number of turns of the SQUID 
superconducting pickup loop around the sample, A is 
the area of the sample. The circuit is schematically 
shown in Fig. [T] Suppose the magnetic field through 
the pickup loop varies sinusoidally as a function of time, 
with frequency cu: B = B^e^'^'^^ . The electromotive 
force around the circuit is then given by Faraday's law: 
£ — iujNABQe~^'^* = iojN^, where we define the mag- 
netic flux $ — BA. The current in the circuit is given by 
/ = 8./{Zin + Zp), where the impedances of the SQUID 
input coil and the pickup loop are in the denominator. 
Expressing these inductive impedances as Z = icuL re- 
sults in / = A^$/(Lin + ip^^), where Li, 



1.5 /iH is the 

input coil self-inductance and Lp^"^ is the self-inductance 
of the N-turn pickup loop. The flux that is coupled to 
the SQUID is given by $SQ = Mn-f, where Mi„ « 10 nH 
is the SQUID-input coil mutual inductance. Thus we 
flnally have: 



SQ 



Li, 



(Al) 



Having a several-turn pickup loop improves the flux cou- 
pling to the SQUID, but only until the pickup coil in- 



ductance reaches L- 



{N) 
P 



If the number of turns is 



increased further, so that Lp^^ ^ Lin, the coupling dete- 
riorates, since Lp^"* scales as iV^. The optimal number of 
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turns can be calculated from Eq. (|Aip , given the dimen- 
sions of the sample and the superconducting wire, and 
usin g a n expression for the inductance of a multi-turn 

coil (sal. 



APPENDIX B: DETUNING OF THE 
GRADIOMETER CONDITION DUE TO SAMPLE 
PERMEABILITY 

Previous work 14] employed two superconducting 
pickup loops around the sample, operated in a concentric 
loop planar gradiometer configuration, to reduce mag- 
netic field noise. Unfortunately, the permeability of the 
sample alters the pickup loop inductances, and spoils the 
gradiometer tuning, and reduces the experiment sensitiv- 
ity slightly because the loop inductance is increased. In 
the case of a single loop or radius a around an infinitely 
long cylinder of permeability jj, and radius b, the increase 
in inductance is [sj] 

/•OO 

dL = 2noa^ ^{k)[Ki{ka)]^dk (Bl) 
Jo 



where Iq and Kq are the zero-order modified Bessel func- 
tions, and the approximation holds when /x >> 1. We 
use SI units in this Appendix. The effect of finite length 
is small, provided that length £ > 2b, which is the case for 
the experiment. In the case a = b for a loop directly on 
a cylinder, SL sa Gfiga w 0.1 fill, for a — 1.2 cm, a 20% 
effect. When a — \/2b, SL = 0.02 /iH, a much smaller 
effect. The loss in sensitivity is around 10%, leaving the 
principal problem due to this effect as the detuning of the 
gradiometer condition, which is temperature dependent, 
and of order 10%. Usually, gradiometers are fabricated 
to provide 1% common mode rejection. Thus, we see 
a substantial limitation to the degree of common mode 
rejection possible when there is a permeable sample in 
the central loop. In principle, the loop areas could be 
adjusted to cancel gradients at a particular temperature, 
but introduces difficulties if temperature variation is em- 
ployed to test for systematic effects. 



where is the permeability of free space, Ki{ka) is the 
ffist order modified Bessel function, and 



<l>{k) = 



[fi ^ l)kbIo{kb)h{kb) ^ Io{kb) 
{fi- l)kbKo{kb)Ii{kb) + 1 ~ Koikb) 



(B2) 
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